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Abstract. The metal content of the spiral galaxy M33 is analyzed through spectroscopic 
observations of its emission-line populations, H n regions and Planetary Nebulae (PNe): 
their abundance gradients are identical within the errors. The 2D metallicity map is pre- 
sented, finding an off-center peak, located in the southern arm. A chemical evolution model 
of M33 described by Magrini et al. (2007a) is updated at the light of recent results on the 
Schmidt law and the metallicity gradients. 
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1. Introduction 

M33 (NGC 598) is a spiral galaxy whose close- 
ness (840 kpc, Freedman et al. 1991; optical 
size 53' x53', Holmberg 1958) and inclina- 
tion (i=53 °) allow detailed studies of its stel- 
lar populations and ionized nebulae. Planetary 
Nebulae (PNe) and H n regions represent two 
very different stages in the lifetime of a galaxy, 
and their comparison provides insight on the 
chemical evolution of the host galaxy. PNe are 
indeed the final ejecta of evolved low- and 
intermediate-mass stars with mass between 1 
and 8 M Q , which must have formed between 
3xl0 7 yr and 10 Gyr ago, e.g., Maraston 
(120051) . while H ii regions belong to a very 
young population. 

One of the most debated questions about 
the chemical evolution of galaxies is the metal- 
licity gradient behaviour with time. Chemical 
evolution models predict different temporal be- 
haviors of the metallicity gradient, depending 



on the assumptions one makes on gas inflow 
and outflow rates, and the star and cloud for- 
mation efficiencies. Observations are needed 
to constrain these assumptions, but so far they 
have been insufficient, especially for the older 
populations. The idea behind the observations 
leading to our work is to study the chemical 
and physical properties of a large number of 
PNe and H n regions in M33, using the same 
set of observations, the same data reduction 
and analysis techniques, and identical abun- 
dance determination methods, to avoid all bi- 
ases due to the stellar vs. nebular analysis. The 
aim is to derive abundances of the a-elements 
for as many PNe and H n regions as possible, 
and to study the variation of the metallicity gra- 
dient and the average abundances in M33. 

In this paper, several aspects of the study 
of M33 will be described. In Section|2]the ob- 
servations of a sample of PNe and H n regions 
with MMT will be presented. These data, to- 
gether with a large literature data-set, allow us 
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to recompute the PN and H n region metal- 
licity radial gradients. In Section [3] the spatial 
distribution of the metallicity will be shown. 
The off-center of the metallicity maps obtained 
from H ii regions and from PNe will be dis- 
cussed. In Section l3~2l the chemical evolution 
model of M33 (Magrini et al. 2007a, hereafter 
M07) will be revised, including a star forma- 
tion process, parametrized by a Schmidt law, 
consistent with the observations. The conse- 
quences, in particular the time evolution of the 
metallicity gradient, will be analyzed. Finally, 
our conclusions and summary will be given in 
Sectiong] 

2. The MMT observations 

In semester 2007B we observed ~150 ionized 
nebulae in M33 in multi-object spectroscopic 
mode. We used the Hectospec fiber-fed spec- 
trograph (Fabricant et al. 2005) on the Multi 
Mirror Telescope (MMT), with a 270 mrrr 1 
grating at a dispersion of 1.2 A pixel" 1 . The 
instrument deploys 300 fibers over a 1° diam- 
eter field of view; the fiber diameter is 1.5" 
(6 pc using a distance of 840 kpc to M33). 
We obtained spectra of 102 PNe and 48 HII 
regions with resulting total spectral coverage 
from 3600 to 9100 A. The details on the data 
reduction and analysis are in Magrini et al. 
(2009, hereafter M09). Most of the PNe ob- 
served in M33 belong to its disk, and they 
are non-Type I, implying a population mainly 
composed by PNe with old progenitors, i.e. 
M<3 M , and ages >0.3 Gyr. 

The MMT observations allowed us to ob- 
tain a good determination of the O/H gradi- 
ent (together with Ne/H, S/H, and Ar/H) both 
from H n regions and PNe. In the case of 
H n regions, our cumulative sample includes: 
i) H ii regions by Magrini et al. (I2007bl) which 
comprises their own determinations and all 
previous abundance determinations with avail- 
able t e , recomputed uniformly; ii) the sample 
by Rosolowsky & Simon (2008); iii) the new 
MMT sample. For PNe, we use the sample of 
PNe presented by M09. The resulting gradients 
are the following, for PNe and H n regions re- 
spectively, 12 + log(0/H) = 

-0.031(+0.013) R GC + 8.44(+0.06), (1) 



-0.032(+0.009) R GC + 8.42(+0.04), (2) 

where Rcc is the galactocentric distance in 
kpc. The two gradient are identical within the 
errors, both in the their slopes and central val- 
ues. 

3. The metallicity in M33 and its 
evolution 

The large amount of chemical abundance data 
from H ii regions and PNe collected to date 
in M33 allow us to reconstruct its metallic- 
ity map. In this section, both the 2D and the 
radial distribution are analyzed taking advan- 
tages also of all previous abundance determi- 
nations with measured t e and of our new re- 
sults. 

3.1. The 2D distribution of metals 

The usual way to study the metallicity dis- 
tribution in disk galaxies is to average it az- 
imuthally, assuming that: i) the centre of the 
galaxy coincides with the peak of the metallic- 
ity distribution; ii) at a given radius, the metal- 
licity is the same in each side of the galaxy. 
The large number of metallicity measurements 
in M33, both from H n regions and from PNe, 
allow us to reconstruct not only their radial 
gradient, but also their spatial distribution pro- 
jected onto the disk. In Figure 1, the two- 
dimensional metallicity distributions for M33 
from H ii regions and from PNe are shown. The 
highest metallicity H n regions and PNe are not 
located at the center of the galaxy, but rather 
lie at a radius of 1 -2 kpc, in the southern arm. 
Simon & Rosolowsky (2008) noticed this be- 
havior for H ii regions and suggested as expla- 
nation that the material enriched by the most 
recent generation of star formation in the arm 
has not yet been azimuthally mixed through the 
galaxy. Whilst it might be true for H ii regions, 
it cannot be the reason for the off-center metal- 
licity distribution of PNe since they belong 
to an o lder population. Colin & Athanassoula 
(1981) noticed several evidences of asymme- 
tries in the inner regions of M33, such as the 
distribution of HI atomic gas, of H n regions, 
of high luminosity stars. They proposed a kine- 
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Fig. 1. The metallicity maps: PNe (top) and H n re- 
gions (bottom). The 12 + log(0/H) scale is shown 
on the y axe. The centre of M33 corresponds to the 
0,0 position. 



matical model with a displaced bulge having a 
retrograde motion around the center. Also de- 
tailed analysis of the innermost regions of M33 
by Corbelli & Walterbos (120071) found possible 
asymmetries in the stellar and gas velocity pat- 
tern which might be related to the displacement 
of a small bulge. The off-center metallicity dis- 
tribution might thus be related to the lack of a 
dominant gravitational source in the center of 
this galaxy with a consequent motion at differ- 
ent epochs of the peak of the highest star for- 
mation region around the M33 visual center. 

3.2. The time-evolution of the 
abundance gradient 

M07 built a chemical evolution model of M33, 
called accretion model, able to reproduce its 
main features, including the radial trends at 
present-time of molecular gas, atomic gas, 



stars, SFR, and the time evolution of the metal- 
licity gradient using the available constraints 
at that time. In that model, the disk of M33 
was formed by continuous accretion of pri- 
mordial gas from the intergalactic medium. 
However, new observations have been made 
recently available, rendering necessary a re- 
vision of the model. In particular, the global 
slope of the radial O/H gradient has been con- 
firmed much shallower than retained in the past 
(cf. Rosolowsky & Simon 2008) and its evolu- 
tion much slower (cf. M09). 

The general assumption of multiphase 
chemical evolution models, such as the M33 
one (M07), is to describe the formation and 
disruption of diffuse gas, clouds, and stars, by 
me ans of physical processes, e.g., Ferrini et 
al. dl992l) . In particular, the SF is represented 
with two processes: the interaction of molec- 
ular clouds with the radiation field of massive 
stars and the collisions between two molecular 
clouds. The dominant process is due to cloud 
collisions. In this kind of model the relation- 
ship between the star formation rate (SFR) and 
the surface density of gas (molecular or total) 
is, thus, a by-product of the model, and cannot 
be assumed 'a priori'. 

However, the law relating the surface den- 
sities of SF and cold gas is one of the most 
fundamental laws describing the galaxy be- 
haviour. Virtually the entire range of global star 
formation rates in galaxies can be reproduced 
by a Schmidt power law relation. In the partic- 
ular case of M33, the relation between the SFR, 
measured from the FUV emission, and molec- 
ular gas has a well-defined slope (Verley et al., 
in prep.) corresponding to 



l. SFR - Al} n ? olgas . 



(1) 



A pure cloud-cloud collision process for the 
star formation is not able to reproduce it. For 
this reason, we have taken into account other 
parameterizations of the SF process. In partic- 
ular, a process which is dominated by cloud 
collisions close to the center, while in the in- 
termediate and peripheral regions is propor- 
tional to the fraction of clouds with a 1.5 ex- 
ponent, is able to reproduce the observations. 
It reproduces the higher cloud surface density 
in the inner regions, rendering the conversion 



4 



Magrini: Metals in M33 




R (kpc) 

Fig. 2. The time-evolution of the O/H radial gra- 
dient. Filled symbols (blue circles and magenta 
squares) are the oxygen abundances of H n regions 
and PNe, respectively, averaged in bins 1 kpc wide. 
Model with the Schmidt law: present time (red con- 
tinuous line), 1 Gyr ago (long-short dashed line), 5 
Gyr ago (dot-dashed line). Model with the cloud- 
cloud collision process (M07): present time (green 
dashed line). 

into star less effective, while in the outer re- 
gions it takes into account a more efficacious 
SF. The resulting Schmidt law would have an 
average exponent all over the radial range (or 
the molecular gas surface density range) of 1 . 1 , 
thus consistent with the observations. 

In addition, the introduction of the Schmidt 
law allowed us to better reproduce the O/H gra- 
dient and its evolution (see Figure 2): a gradi- 
ent almost flat both at present-time both at the 
epoch of the formation of the PNe progenitors, 
with a little evolution of its absolute value and 
slope. Note from Figure 2 the steeper gradient 
predicted by the previous model, where the SF 
process was dominated by cloud-cloud colli- 
sion all over the radial range. 

4. Summary and conclusions 

The chemical evolution of M33 is studied by 
means of new spectroscopic observations of 
PNe and H n regions. Their 2D metallicity 
maps have been found both off-centered, with 



a peak in the southern arm, at 1-2 kpc from 
the center. This might be related to the absence 
of a dominant gravitational source in the cen- 
ter. The slow evolution of the metallicity gra- 
dient from the present time to the birth of the 
PNe progenitors is explained with an accretion 
model where the SF process is driven by the 
Schmidt law. 
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